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We present results of coupled-channels calculations which are performed to understand 
the role of vibrational excitations of the daughter nucleus in proton emission decay of spher- 
ed , ical nuclei. It is shown that the experimental decay rates from the 2d 3 / 2 , 3si/2 and 
(— i ■ states in proton emitting nuclei 160,161 Re are reproduced simultaneously when the quadrupole 
phonon excitation of the core nuclei is taken into account. We apply the formalism also to 
the 145 Tm nucleus and show that good agreement is achieved with the recently measured 
branching ratio for the fine structure. 

t-H \ 

§1. Introduction 

>. 

Although nuclei beyond the proton drip line are unstable against proton emis- 
sion, their lifetime is sufficiently long to study the spectroscopic properties, because 
of the fact that the proton has to penetrate the Coulomb barrier. Thanks to the 
recent experimental developments of production and detection methods, a number 
ground-state as well as isomeric proton emitters have recently been discovered. EP'13 
These experimental data have indeed shown that the proton radioactivity provides 
a powerful tool to study the structure of proton rich nuclei. 

In the field of heavy-ion reactions, the coupled-channels formalism has been 
a standard framework to extract nuclear structure information from experimental 
data. This formalism has been successfully applied also to deformed proton emit- 
ters, i^'i^'i) For instance, the (52 deformation parameter of deformed proton emitters 
has been extracted from the experimental decay widths a/o branching ratios. 

In the previous analyses for the proton emitters in the region 68 < Z < 82, where 
proton emissions from the l/in/2> 3 s i/2i an d 2d 3 / 2 orbitals have been observed, the 
emitters have been treated as spherical nuclei without any structure. & It has been 
pointed out that such calculations systematically underestimate the measured decay 
half- lives for proton emissions from the 2^3/2 state, while the same model works well 
for emissions from the l/in/ 2 and 3si/2 states in both odd-Z even-N and odd-Z odd- 
N nuclei. For the 151 Lu nucleus, this discrepancy was attributed to the effects of 
oblate deformation of the core nucleus^ 50 Yb, which alter both the decay dynamics 
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and the spectroscopic factor. e])>EJ)'II2P This fact has motivated us to investigate 
the effects of vibrational excitations on the proton decay of spherical nuclei. 0^ '111* 
Similar effectshave been well recognized in the field of heavy-ion subbarrier fusion 
if,® 



reactions. 



In this contribution, we solve the coupled-channels equations for spherical proton 
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emitters in order to study whether the effects of vibrational excitations of the daugh- 
ter nucleus consistently account for the measured decay half-life from the lh n / 2 , 
3si/2, and 2d 3 / 2 states. In Sect. 2, we briefly review the coupled-channels frame- 



work for proton emissions. We use the Gre 



unction approach to compute the 



decay rate of resonance states. Q)'EP'lla)'E3)'t3)'llHP We discuss the method both for 
the one dimensional and for the coupled-channels problems. In Sect. 3.1., we com- 
pare our results with the experimental data. We particularly study proton emi ssio ns 



from the lh n / 2 and 3sy 2 states of 161 ReEl and from the 2d 3 / 2 state of 160 Re.EJ A 
clear vibrational structure in the spectrum of the core nucleus 160 W for the proton 
emitter 161 Re has recently been identified up to the state.© We shall assume 
that the vibrational properties are identical between the core nuclei 160 W and 159 W, 
and study the sensitivity of the calculated decay rate to the dynamical deformation 
parameter for the vibration. We will also discuss the dependence on the multipolar- 
ity of the phonon state. In Sect. 3.2., we apply the formalism to the proton decay 
of the 145 Tm nucleus, where the fine structure was recently observed. c$ We then 
summarize the paper in Sect. 4. 

§2. Green's function method for proton radioactivities 

2.1. One dimensional case 

In this section, we review the method which we use below to calculate the decay 
rate for proton emissions. Let us first consider the no coupling case. Extension to 
the coupled-channels problems will be given in the next subsection. 

We consider the following Schrodinger equation for a resonance state: 



h 2 d 2 l(l + i)h 2 . . f„%„ 

+ V(r) -{E + -r 



2/idr 2 2/xr 2 V 2 



u(r) = 0, (2-1) 



where \i is the reduced mass for the relative motion between the valence proton and 
the daughter nucleus, and V(r) is a potential which consists of the nuclear, the spin- 
orbit, and the Coulomb parts. We describe the resonance state using the Gamow 
wave function. Accordingly, u(r) satisfies the boundary condition that it is regular 
at the origin and is proportional to the outgoing wave asymptotically, namely 

u(r) = r l+1 (r -> 0), (2-2) 

-» N{Gi{kr) + iF l {kr)) (r -»■ oo), (2-3) 

where k is the wave number, and F\ and G\ is the regular and the irregular Coulomb 
wave functions, respectively. As is well known, these boundary conditions cannot be 
satisfied with a real energy E, and an imaginary part of the energy is included in 
Eq. fl2-l|) . The latter is nothing but the decay width of the resonance state. This 
quantity is calculated also as a ratio between the asymptotic outgoing flux and the 
normalization inside the barrier, EP 



in — Wf \u(r)\ 2 dr. (2.1) 
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where r< is the outermost turning point. 

An alternative approach to calculate the decay^radih En iias_been proposed, 
which is based on the Green's function formalism. cP' 13* >ll9)'EJ)'II^'[L£P This method 
was first developed for a decays by Kadmensky and his collaborators © and was 
recently applied to the coupled-channels problem for deformed proton emitters by 
Esbensen and Davids. o> 

In this method, one first sets i~b to be in Eq. ( p-l| ) and looks for a standing 
wave solution which has an asymptotic form of 



</>(r) ->jVGi(kr) 



(r 



oo 



(2-5) 



The coefficient N is determined so that the wave function <f> is normalized to one 
inside the barrier. Note that (j) corresponds to the solution where the nuclear phase 
shift is equal to tt/2. This procedure should be a good approximation since i~o is 
extremely smaller than the real part of the energy E (typically i~o is between 1CP 18 
and 10" 22 MeV while E is order of MeV) . % The next step is to apply the Gell-Mann- 
Goldberger transformation^^ to the total Schrodinger equation [T — (E +ir$ /2)] l P = 
—V&, where T is the kinetic energy operator. Adding the point Coulomb potential 
Vc(r) = Zf)e 2 /r, Zd being the charge number of the daughter nucleus, to both the 
left and the right hand sides, this equation is transformed to 



T + 



Z D e 2 



'Z D e 2 



V(r) V. 



The formal solution of this equation reads 



Z D e 2 



T + 



E+i^ 



V(r) 



(2-6) 



(2-7) 



Next two approximations are made. Firstly one replaces on the right hand side of 
Eq. ( |2-7|) with the standing wave solution 0, whose radial part satisfies the boundary 
condition given by Eq. (|2-5| ). Secondly, one regards i~o/2 as an infinitesimal number r\ 
which provides the outgoing wave boundary condition for the single particle Green's 
function (T + Vc — E — ir/)^ 1 . As long as i~b is extremely small compared with E, 
its actual value is not important at this stage. As is well known, the single particle 
Green's function is expressed in terms of the regular and the outgoing wave functions 



as 



19) 



1 



T + Vc-E-ii] 



1^^— V~ yjlm(r>)yjim(r<), (2-8) 



Lm 



where 0/ = Gi + %F\ is the outgoing Coulomb wave. One then obtains 





kh 2 



T + Vc-E-ir) 



yjlm(r) / dr' Ftikr') 



Z D e 2 



V(r') 0, (2-9) 



Z D e 2 



V(r'))ct>(r'), (2-10) 
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as r — > oo. Prom this equation, the normalization factor J\f in Eq. (2-3) is found to 
be 

M = wC dr ' Fl[kr ' ] i^v~ ~ v{r,) ) 0(r ' } - (2 ' n) 



The decay rate is then obtained by Eq. (2-4), by replacing u{r) with (f>(r) in the 
denominator. 

2.2. Coupled- channels problem, 

Let us now discuss the channel coupling effects on proton radioactivities of spher- 
ical nuclei. In order to take into account the vibrational core excitations, we consider 
the following Hamiltonian: 

— V 2 
2fi 



H = -— V 2 + V(r) + V coup (r, a) + H vih , (2-12) 



where a is the coordinate for the vibrational phonon of the daughter nucleus. It is 
related to the dynamical deformation parameter as = /^x+i ( Q V + ( — ) Ma A^)> 
where A is the multipolarity of the vibrational mode and (a^) is the creation 

(annihilation) operator of the phonon. H vl ^ = hoj a \fj, a \^ is the Hamiltonian 
for the vibrational phonon. In this paper, for simplicity, we do not consider the 
possibility of multi-phonon excitations, but include only excitations to the single 
phonon state. The coupling Hamiltonian V coup (r,a) consists of three terms, i.e., 

V CO up(r,a) = Vc<Zp(r,a) + Vcon P (r,a) + Vcou P (r,a). The nuclear term reads 

where the dot denotes a scalar product. We have assumed that the nuclear potential 
is given by the Woods-Saxon form. VXr(r) is the nuclear component of the bare 
potential V[r). This term is subtracted in order to avoid the double counting. Notice 
that we do not expand the nuclear coupling potentiaL hiijt include the couplings to 
the all orders with respect to the phonon operator. u>'c3> The matrix elements of 
this coupling Hamiltonian are evaluated using a matrix algebra, as in Ref. ^). As 

for the Coulomb Vmup as well as the spin-orbit Vcoup terms, the effects of higher 
order couplings are expected to be small, c3) and we retain only the linear term. The 
Coulomb term thus reads 

V^ ) (r,a) = ^^^ Tl (^) X a x .Y x {r) (for r > R c ) (2-14) 



r 

o2 i / „ \ A 



3Z D e z 1 / r 



R c 2A + 1 \R. 



a x -Y x (r) (for r < R c ), (2-15) 



where R c is the charge radius. For the spin-orbit interaction, we express it in the so 
called Thomas form,o''E§' 

V«%(r, a )=V so - r f r l-<T + i V so R so £ { (V|; Y XfM (r )) • (V x <x)j , (2-16) 
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where f(r) = 1/[1 + exp((r — R so )/a so )]. The last term in Eq. ,( 12/16 ) can be decom- 
posed to a sum of angular momentum tensors using a formula 



27F 



In order to solve the coupled-channels equations, we expand the total wave 
function as 

,(?) 



U i"pj p nl( r ) 

7 



# jm (r,a) = EE I(Wp^'™>> (2-18) 



where 



{r,a\(l p j p nl)jm) = ^2 (j p m p Im I \jm)y jp i pmp (r)ip nImi (a), (2-19) 

99 being the vibrational wave function. The coupled-channels equations are obtained 
by projecting out the intrinsic state \(l p j p nl)jm) from the total Schrodinger equa- 
tion, and read 



h 2 d 2 UL + l)h 2 , » / i . 



u 



Ipjp 



nl( r ) 



= ~ E E((^" J )j' m l^oup|(^> / lOim}wg, ri , J ,(r)(2-20) 

To solve these equations, one needs to compute the coupling matrix elements of 
the operators a\ ■ T\, where T\ M is either or_jY^±i (— i V x <x)]( A/x ). These are 
expressed in terms of the Wigner's 6-j symbol asB3) 



((l'j' p n'l')jm\a x ■ T x \(l p j p nl)jm) = 



3 I' f P 
A j p I 



2-21) 



^(yj'pl'p\\Tx\\yjplp)(<Pn'I'\\ax\\^nl)- ( 

For transitions between the ground and the one phonon states which we consider in 
this paper, the reduced matrix element (^'n'/'II^All^'n/ 

} is given by (3\. The reduced 
matrix elements for the operators T\ are found in Ref. 28). 

The Green's function method discussed in the previous subsection can provide 
a convenient way to calculate the decay width particularly for the coupled-channels 
problems. & The direct method seeks the Gamow wave functions where uf^ nI {r) 

have the asymptotic form of J^i p j pn j(Gi p (k n jr) + i Fi p (k n jr)) at r — > 00 for all the 

channels, where k n i = \J2/j(E — nTwj)/Ti 2 is the channel wave number. This method, 
however, requires to solve the coupled-channels equations in the complex energy 
plane and out to large distances, which is quite time consuming and also may be 
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difficult to obtain accurate solutions. In the Green's function method, in contrast, the 
coupled-channels equations are solved in the real energy plane and the solutions are 
matched to the irregular Coulomb wave functions G[ at a relatively small distance 
r matc h, which is outside the range of nuclear couplings. From the solution of the 
coupled-channels equations <F^(r, a) thus obtained, the outgoing wave function for 
the resonance Gamow state is generated using the Coulomb propagator as&'lll> 



&j m (r, a) = — J dr'da (ra 



1 



-ffcoul + #vib — E — IT] 



r'ol 



x (V(r') + WV) - ^f] *%n(r',a'), (2-22) 



where H cou \ = —h 2 S/ 2 /2Li + Z^e 2 /r is the Hamiltonian for the point Coulomb field. 

r(j) 

IpjpTlI 



r(i) 

As in the previous subsection, the asymptotic normalization factors M, ■ T are found 



to beB) 

2// 




V(r) + V cmp (r,a) 



r 



( 2 - 23 ) 



h 2 k n i Jo 



drF lp (k nI r) { [V(r) - ) u lpjpnI (r) 



+ E E((w P ^)i^i^ou P i(/;j> , /Oj^)^gv„'/'(0 (2-24) 

In this way, the effects of the long range Coulomb couplings outside the matching 
radius r matc h are treated perturbatively. The partial decay width is then calculated 
as ripjpni = F^knilNipjpnil 2 / provided that the wave function is normalized 
inside the outer turning point. 

§3. Comparison with experimental data 

We now solve the coupled-channels equations and compare the results with the 
experimental data. One of the most important experimental observables is the decay 
half-life. In the present framework, it is calculated as 

T 1/2 = ^ln2, (3-1) 

where Sj is the spectroscopic factor for the resonance state, which takes into account 
correlations beyond the single-particle picture. If one assumes that the ground state 
of an odd-Z nucleus is a one-quasiparticle state, the spectroscopic factor Sj is iden- 
tical to the unoccupation probability for this state and is given by Sj = u 2 in the 
BCS approximation. %> 
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In all the calculations in this sec- 
tion, we use the real part of the 
Becchetti-Greenless optical model po- 
tential forthe proton-daughter nucleus 
potential. c3 The potential depth was 
adjusted so as to reproduce the exper- 
imental proton decay Q value for each 
value of the dynamical deformation pa- 
rameter (3 and the excitation energy Two 
of the vibrational phonon excitations in 
the daughter nucleus. Following Ref. 
0), we assume that the depth of the 
spin-orbit potential is related to that of 
the central potential by V so = — 0.21-0 • 
The charge radius R c is assumed to be 
the same as R$ in the nuclear poten- 
tial. The spectroscopic factor Sj in 
Eq. ( |3-1| ) is taken from Ref. |^). This 
was evaluated in the BCS approxima- 
tion to a monopole pairing Hamiltonian 
for single-particle levels obtained with a 
spherical Woods-Saxon potential. 

3.1. 160,161 i?e nuclei: reconciliation of 
the 0I3J2 puzzle 

In this subsection, we solve the 
coupled-channels equations for the res- 
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Fig. 1. The decay half- lives for proton emis- 
sion from the 2d 3 / 2 state of 160 Re and the 
3sx/2 an d lftn/2 states of 161 Re as a func- 
tion of the dynamical deformation param- 
eter P2 of the quadrupole vibrational ex- 
citation of the daughter nuclei. The vi- 
brational excitation energy is hu>2 = 0.61 
MeV. The experimental data are taken 
from Refs. 20) and 21) and are denoted 
by the dashed lines. The arrows indicate 
the results in the no coupling limit. 



as well as the 2d 3/ / 2 state in 160 Re. The 
first 2 + state in the core nucleus 160 W for the 161 Re nucleus was recently observed 
at fiuJ2 =610 keV.EHP Accordingly we first discuss the effects of the quadrupole vibra- 
tional excitations. Figure 1 shows the dependence of the decay half-life for proton 
emissions from the three resonance states on the dynamical deformation parame- 
ter fa °f the quadrupole mode. The experimental values for the decay half-life are 
taken from Refs. |20|),^l|), and are denoted by the dashed lines. The arrows are the 
half-lives in the absence of the vibrational mode. For the decays from the 3si/ 2 and 
1/&1X/2 states, if one takes into account uncertainty in the experimental Q value of the 
proton emission, these calculations for the decay half-life in the no coupling limit are 
within the experimental error bars.U In clear contrast, the calculation for the 2g?3/ 2 
state is still off from the experimental data by_about 30% even when uncertainty of 
the decay Q value is taken into consideration.™ 

The results of the coupled-channels calculations are shown by the solid lines in 
the figure. We treat the extra neutron as a spectator for the odd-odd proton emitter 
160 Re. One notices that the channel coupling effects significantly enhance the decay 
half-life for the 2d 3 / 2 state, while the effects are more marginal for the 3si/ 2 and 
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l/iXi/2 states (notice the difference of the 
scale of the vertical axis in the figure). 
Our results show that all of the three 
half-lives can be reproduced simultane- 
ously if the dynamical deformation pa- 
rameter is around fa ~ 0.16. 

Why is only the d 3 / 2 state affected 
considerably by the channel coupling ef- 
fects while the others are not? A pos- 
sible explanation is as follows. The 
quadrupole interaction mixes the s and 
d states. For the S]/ 2 state, a small ad- 
mixture of d states does not affect much, 
since the decay from the |si/2<8>0 + ) state 
dominates after all. On the other hand, 
for the dz/2 state, a small admixture of s 
state can influence the decay dynamics 
significantly depending on the excitation energy: there is no centrifugal barrier for 
the \s\/2 <8> 2 + ) state, and thus the decay from this component can compete with the 
decay from the dominant component ^3/2 ® + ) state even though the energy for 
the relative motion is smaller. For the /in/2 state, the quadrupole interaction mixes 
it with the I/7/2 ® 2 + ) component. The decay from the latter may compete with the 
decay from the main component due to the smaller centrifugal potential. Compared 
with the c?3/2 state, however, the effective proton energy in the incident channel de- 
fined as Ep" = E p — V cen t(r) — V(r) is smaller, which makes the channel coupling 
effects weaker. Therefore, the core excitations do not influence significantly the pro- 
ton decay from the foii/2 state unless the coupling is very large a/o the excitation 
energy is very small. 

Let us now discuss the dependence of the decay rate on the multipolarity of 
the phonon mode. Figure 2 shows the effects of octupole phonon excitations on the 
decay half-life from the 2d 3 / 2 state of 160 Re. As an illustrative example, we take 
^3=0.2 MeV, but results are qualitatively the same for different values of T10J3. For 
the octupole vibration, the enhancement of the half-life is too small to account for 
the observed discrepancy between the experimental decay half-life and the prediction 
of the potential model with no coupling. As was noted before, IHP the proton decay is 
very sensitive to the angular momentum of the proton state. The same seems to be 
true for the multipolarity of the collective excitation of the daughter nucleus. 

3.2. 145 Tm nucleus: fine structure 

We next discuss the fine structure in proton emission from the 145 Tm nucleus, 
which was recently observed at the Oak Ridge National Laboratory.© In the exper- 
iment, two clear peaks were seen in the energy spectrum of the emitted proton, one 
corresponding to the transition from the l/in/2 state in 145 Tm to the ground state 
of the daughter nucleus 144 Er and the other to the first excited state at 0.326 MeV. 
The measured ground state decay half-life was 3.1±0.3/is, and the branching ratio 




Octupole 
0)3= 0.2 MeV 



400.' 1 1 1 1 1 1 

X).l 0.15 0.2 0.25 

Fig. 2. The decay half-life for proton emission 
from the 2(1^/2 state of 160 Re as a function 
of the dynamical deformation parameter of 
the octupole phonon excitation of the core 
nucleus. The excitation energy is set to be 
0.2 MeV. 
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was Davids and Esbensen 

have made a detailed analysis for 
the fine structure of this decay mode 
using the particle-vibration coupling 
model. Hif They obtained good agree- 
ment with the data for both the half- 
life and the branching ratio. Figure 3 
shows results with our model Hamilto- 
nian, which is slightly different from the 
one used by Davids and Esbensen (e.g., 
we include the coupling Hamiltonian for 
the spin-orbit interaction as well). The 
experimental data are denoted by the 
dot-dashed lines. The three theoreti- 
cal curves are for three different pro- 
ton energies, taking its experimental un- 
certaintity into consideration. We find 
that the ground state half-life and the 
branching ratio are simultaneously re- 
produced if the dynamical deformation 
parameter 02 is larger than 0.22, al- 
though more studies on the dependence 
of the observables on the strength of 
the spin-orbit interaction V^ would be 
needed in order to extract a more con- 
clusive value of 02- Our results indicate 
that the |/ 7 / 2 <8> 2 + ) state admixes 1.08% 
® + ) state at 02 = 0.18. This is co 
given in Ref j23|) . 




2 1 , 1 , 1 , 1 l 




P 2 

Fig. 3. Sensitivity of the decay half-life (the 
upper panel) and the branching ratio (the 
lower panel) for the proton emission of 
145 Tm to the dynamical deformation pa- 
rameter 02 of the quadrupole phonon exci- 
tation of the daughter nucleus 144 Er. 

in probability with the main component 
isistent with the experimental suggestion 



§4. Summary 



We have solved the coupled-channels equations to take into account the effects 
of the vibrational excitations of the spherical daughter nucleus during the proton 
emission. Applying the formalism to the the resonant l/in/2 and 3s\/2 states in 
161 Re and the 2^/2 state in 160 Re, we have found that the experimental data for 
the decay half-lives for these three states can be reproduced simultaneously if the 
quadrupole phonon excitation with 02 ~ 0.16 is considered. This removes the dis- 
crepancy observed before between the experimental data and the prediction of the 
optical potential model calculation for the decaying 2d 3 / 2 state in this mass re- 
gion without destroying the agreement for the l/in/2 anc ^ ^ s i/2 states. A similar 
calculation with the octupole phonon was not satisfactory. We also discussed the 
fine structure in proton emission of 145 Tm nucleus. The experimental data for the 
half-life as well as the branching ratio were simultaneously reproduced within our 
model, with a reasonable admixture of the excited state components in the total 
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wave function. 

The previous coupled-channels analyses were restricted to deformed proton emit- 
ters. Our studies indicate that the channel coupling effects are significant also for 
spherical emitters. Combining with the previous results, our analyses show that all 
the observed proton emitters can now be described in a simple single particle picture 
provided that the collective motion of the daughter nucleus is taken into account. 
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